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The reaction ^''^Pr(Q,n)^''''Pm was investigated between Ea — 11 MeV and 15 MeV with the 
activation method using the 77 coincidence method with a segmented clover-type high-purity Ger- 
manium (HPGe) detector. Measurements with four other HPGe detectors were additionally made. 
The comparison proves that the 77 coincidence method is an excellent tool to investigate cross 
sections down to the microbarn range. The (a,n) reaction at low energy is especially suited to 
test a-|-nucleus optical-model potentials for application in the astrophysical p process. The experi- 
mentally determined cross sections were compared to Hauser-Feshbach statistical model calculations 
using different optical potentials and generally an unsatisfactory reproduction of the data was found. 
A local potential was constructed to improve the description of the data. The consequences of ap- 
plying the same potential to calculate astrophysical (7,0) rates for ^*^Pm and ^^*Gd were explored. 
In summary, the data and further results underline the problems in global predictions of a-|-nucleus 
optical potentials at astrophysically relevant energies. 

PACS numbers: 25.55.-e, 26.30.-k, 27.60.-|-j 



I. INTRODUCTION 

According to the current understanding most nuclei 
heavier than iron are synthesized by neutron-capture re- 
actions via the s and r processes [I-Q . About 35 proton- 
rich nuclei in this mass region, however, are bypassed by 
these processes. These nuclei are called p nuclei. The 
origin of the p nuclei is not completely understood and 
contributions of different independent processes for their 
production are under discussion. The original suggestion 
to produce proton-rich nuclides in a p process, (i.e., by 
proton capture reactions in the H-rich envelope of type 
II supernovae [1]) was later shown to be unfeasible 0, Sl- 
it was found, however, that the 0/Ne layers of a mas- 
sive star are sufficiently heated during the passage of the 
explosive shock wave of a core-collapse-induced super- 
nova to allow partial photodisintegration of the s- and 
r-process nuclei previously present in the stellar plasma 
[9l-[ll|. The photodisintegration reactions produce p nu- 
clei through sequences of (7,n), (7,p), and (7,0;) reactions 
at plasma temperatures in the range of 2 < T < 3 GK. 
This so-called 7 process is the currently accepted nucle- 
osynthesis mechanism to explain the majority of p nuclei. 
Among the few exceptions are ^"^^La and ^^°™Ta, both 
of which have very low abundances in the solar system. 
They could be produced through neutrino reactions with 
neutrinos emitted by the nascent neutron star emerging 
from the core collapse (the v process) Some p nu- 

clei (i64Er and ^""^GA [11, "^in and "'^Sn may also 
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receive stronger contributions from the s and r process 
than previously estimated and may not require a large 
7-process production. 

In addition to the above exceptions, two mass regions 
have remained problematic in explaining the production 
of p nuclei by the 7 process in core-collapse supernovae: 
the lightest p nuclei with mass numbers A < 100 and 
those in an intermediate region at 150 < A < 165 are 
under produced [1, [l^, [HI, [H, [3 ■ I* ^^^^y un- 

derstood to which extent the deficiencies found are due 
to the astrophysical modeling or the nuclear physics in- 
put (reaction rates) [itI]. For example, recent simulations 
[l8| found that light p nuclei are produced in sufficient 
amounts in a 7 process occurring in the thermonuclear 
explosion of a white dwarf (type la supernova), in con- 
trast to earlier simulations finding no such production 
[19I - I22I . Regardless of the site, such calculations also im- 
plement a 7 process and require a sound determination of 
the astrophysical reaction rates appearing from nuclear 
physics. 

In total, the 7 process involves an extensive reaction 
network consisting of about twenty thousand reactions 
on approximately two thousand nuclei. Most of these 
nuclei are unstable and therefore are not directly accessi- 
ble for experiments. Due to the astrophysically relevant 
low interaction energies, data is also scarce in the rele- 
vant energy range for stable nuclei. In consequence, reac- 
tion network calculations for the 7 process are based al- 
most completely on theoretically predicted reaction rates 
stemming from Hauser-Feshbach statistical model calcu- 
lations [ill, [iSl • The accuracy of the predictions mainly 
depends on the adopted nuclear models for optical-model 
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potentials, 7-strength functions, and nuclear-level densi- 
ties. 

The low-energy a+nucleus optical potential deter- 
mines deflections in the photodisintegration path at in- 
termediate and high masses and therefore impacts the 
calculated p abundances also in the problematic region 
of 150 < ^ < 165 [H [ii]. The few experimental stud- 
ies available close to the astrophysical energy region in 
this mass range revealed a systematic overprediction of 
the (a, 7) cross sections with the widely used potential of 
[27| (see, e.g., (28l - l33 | and references therein). The pre- 
dictions are mostly factors of 2 to 3 above the data, with 
the exception of [35| where the measured ^''^Sm(a,7) S 
factors were found to be lower by more than an order 
of magnitude than the standard prediction at astrophys- 
ical energies. So far, the data at low energy are still too 
scarce to allow the construction of a global optical poten- 
tial suited to predict astrophysical reaction rates further 
off stability, as required by the 7 process. 

In the present work, the reaction ^^^Pr(a,n)^^^Pm 
was studied with the activation technique (for details, 
see [S^l) at the cyclotron of the Physikalisch-Technische 
Bundesanstalt (PTB) in Braunschweig, Germany [s^ to 
improve the experimental situation on the a-|-nucleus 
optical-model potential. The case of ^^^Pr-t-a is espe- 
cially interesting because ^*^Pr is close to -'^^^Sm and is 
also neutron-magic. It is important to see whether sim- 
ilar problems in the prediction of cross sections and as- 
trophysical S factors arise for reactions on this nucleus 
as those found for ^'*^Sm. Sensitivity studies show that, 
in the case of ^^^Pr, the (a,n) reaction is better suited 
to improve the a-l-nucleus optical-model potential com- 
pared to the (a, 7) reaction (see Sec. [iTl- The (a,n) re- 
action was investigated at laboratory a energies between 
11.0 MeV and 15.0 MeV. 

The sensitivities of the predicted reaction cross sec- 
tions to different nuclear input are discussed in Sec. |lll 
After presenting the experimental method in Sec. Illll the 
data analysis is explained in Sec. |IVl The deduced cross 
sections are compared to theoretical predictions from 
Hauser-Feshbach statistical model calculations in Sec. El 



II. SENSITIVITY STUDIES 

As described in the introduction, photodisintegration 
reactions such as ^^^Pm(7,a)^^^Pr play a major role in 
the 7 process. A disadvantage common to all direct pho- 
todisintegration experiments is the fact that these mea- 
surements, in general, can only account for transitions 
stemming from the ground state of the target nucleus 
whereas the high temperatures in an explosive astrophys- 
ical environment lead to a significant thermal population 
of excited levels and reactions proceed on nuclei in ex- 
cited states as well. The modification of the astrophysical 
reaction rates due to the thermal population is strong, 
especially at temperatures typical for explosive nucle- 
osynthesis. Stellar photodisintegration rates can differ 
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FIG. 1. Sensitivities s of the ^*^Pr(Q,7) laboratory cross sec- 
tions when varying neutron, proton, a, and 7 widths sepa- 
rately by a factor of 2. The astrophysically relevant energy 
range for T = 3 GK is marked by the shaded area. Within 
this energy window the cross sections are almost exclusively 
sensitive to the a width, whereas at energies measurable with 
the activation method within a reasonable time (well above 
10 MeV), the cross section prediction is sensitive to additional 
nuclear parameters as well. 

by several orders of magnitude from laboratory ground- 
state rates It can be shown, however, that the stellar 
effects are smaller in reactions with positive reaction Q 
value as compared to their counterparts with negative Q 
value [13, • Furthermore, it is a commonly adopted 
method when using fits of reaction rates in stellar mod- 
els to only include the reaction direction with positive Q 
value and to compute its reverse rate by application of 
detailed balance [H, [3^-[4l| . This ensures the fit accu- 
racy and numerical stability in the solution of the coupled 
differential equations of the reaction network. 

Therefore, the preferred experimental strategy for as- 
trophysical implementation would always be to measure 
reactions with positive Q values. In this manner, mea- 
surements can study cross sections as close as possible 
to those required for the stellar reaction rates. There 
are comparatively few exceptions to this rule but they 
always enforce the importance of measuring captures in- 
stead of photodisintegrations [H,!!^. For example, both 
^^^Pr(a,7) and ^''^Sm(a,7) have negative Q values but 
the stellar effects in their rates are smaller than in their 
reverse rates. 

Furthermore, it is very important to measure in the as- 
trophysically relevant energy range or at least as close as 
possible to those energies. This ensures that dependen- 
cies of the predictions on nuclear properties are similar 
to those appearing in the stellar rates because these sen- 
sitivities strongly vary with energy. 

The above considerations are also essential in the re- 
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FIG. 2. Sensitivities s of the ^*^Pr(ci,n) laboratory cross sec- 
tions when varying neutron, proton, a, and 7 widths sepa- 
rately by a factor of 2. 



action studied in this work. Although the ^'''^Pm(7,Q;) 
reaction would act at the high plasma temperature in 
the 7 process, it is more advantageous to experimentally 
study the inverse reaction ^'*^Pr(a,7)^'*^Pm. The sensi- 
tivities of its laboratory cross sections to a change of the 
different input parameters are shown in Fig. [T] 

A detailed discussion of the usefulness of the sensitivity 
factor can be found in Ref. js^]. The sensitivity factor s 
describes the change in the cross section when one of the 
nuclear physics input parameters is changed by a factor 
of / = r'/r, where F is the width (or transmission coef- 
ficient) before a variation and F' the modified width. A 
sensitivity s = 1 means that the cross section is changed 
by the same factor as the input parameter, while s = 
signifies that the cross section is not changed at all; that 
is, the input parameter has no infiuence on the predicted 
cross section. Thus, if a width is changed by a factor 
/ = F'/F, the cross section will change by 



the a-|-nucleus optical-model potential. At energies mea- 
surable with the activation method within a reasonable 
time (well above 10 MeV), the cross-section prediction is 
additionally sensitive to the 7 and neutron widths. Thus, 
it is not easily possible to extrapolate the experimental 
data towards the astrophysically relevant region without 
further assumptions. Furthermore, it would be difficult 
or impossible to disentangle the impact of the different 
sensitivities if a discrepancy between experiment and pre- 
dicted cross sections were found. 

For this reason, we chose to perform an experiment 
using the reaction -'^^-'^Pr(a,n)^^''Pm instead. Because of 
its larger cross section, a shorter half-life of the reaction 
product ^^*Pm, and stronger 7 intensities in ^''^Nd [com- 
pared to the (a, 7) reaction and its reaction and decay 
products] , it is possible to measure this reaction down to 
lower energies. The sensitivity of the {a,n) cross section 
with variation factors / = 2 is illustrated in Fig. Ex- 
cept close to the (Q;,n) threshold, the cross sections are 
only sensitive to the a width. Thus, a comparison of the 
measured cross sections with calculated ones allows us 
to study the a-l-nucleus optical-model potential even at 
energies where this is not possible by a capture. Unfor- 
tunately, the (a,n) threshold is located above the upper 
edge of the astrophysical energy window of the capture 
reaction and therefore the measurements cannot be ex- 
tended into this region. Nevertheless, the a-f nucleus op- 
tical potential can be investigated at lower energies than 
previously available for this target. A previous measure- 
ment determined (a,n) cross sections with the activation 
method only at higher energies, between 15 and 45 MeV 
[4^. In addition, their lowest energy data point at 15.71 
MeV carries a very large uncertainty on the energy and 
cannot constrain the predictions. Within this uncertain- 
ties it agrees with our data shown in Sec. [V] 

It should be mentioned that a-l-nucleus potentials can 
also be studied in elastic scattering experiments [46|. 
However, these experiments sometimes suffer from am- 
biguous solutions 47 1 . 



III. EXPERIMENT 
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In Fig. [T] the neutron, proton, a, and 7 widths were 
each varied independently by a factor of / = 2. The 
energy window relevant for the calculation of the reac- 
tion rate is located between 7.2 MeV and 10.3 MeV for 
T = 3 GK [i^l. In this energy region the cross section 
prediction is almost exclusively sensitive to the variation 
in the a width. This width, in turn, is determined by 



All activations were carried out at the cyclotron of the 
Physikalisch-Technische Bundesanstalt (PTE) in Braun- 
schweig [36| . Thin samples were irradiated for several 
hours and the 7 rays emitted after the decay of the ra- 
dioactive reaction products were detected with a clover- 
type high purity Germanium (HPGe) detector at the In- 
stitut fiir Kernphysik in Cologne using the 77 coincidence 
method. 



A. Preparation and characterization of targets 

Metallic praseodymium was evaporated in vacuum 
onto 1-mm-thick aluminum disks with a diameter of 35 
mm. These backings were thick enough to stop the a 



4 



suppression aperture 
I voltage; Us= -300V current 



s 



s 



■3 



1.8 
1.6 
1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



in 



• Experimental data from RBS measurement 



irradiated area 
_] I i_ 



5 10 15 20 

Position on Target [mm] 



FIG. 3. Homogeneity of the target. The areal particle den- 
sity of this target was measured at nine positions distributed 
over the whole target. Two measurements were performed 
on the Al backing; four on the part onto which Pr was evap- 
orated but which was not irradiated with q particles; three 
measurements were performed on the irradiated part of the 
praseodymium. These regions are marked in white, light gray, 
and dark gray, respectively. The depicted target was irradi- 
ated in the activation experiment with a particles of Ea = 11 
MeV. The other targets used showed similar homogeneities. 



beam completely to ensure a reliable charge collection. 
The Pr samples were characterized by Rutherford back 
scattering (RBS). The RBS measurements were carried 
out with singly charged *He ions with an energy of 2 
MeV ± 1 kcV at the RUBION Dynamitron-tandem ac- 
celerator laboratory at the Ruhr University of Bochum. 
The targets were mounted on a movable sample holder, 
which serves as a Faraday cup and were irradiated with 
a beam current of 14 uA. Secondary electrons were sup- 
pressed by an appropriate voltage (—300 V) to ensure 
a reliable charge collection using a calibrated integrator. 
The silicon detector was mounted at a distance of 35 mm 
from the center of the sample. Its resolution is 15 keV 
for 2 MeV a particles. The solid angle of the detector 
was determined to be (1.91 ± 0.07) msr. The areal par- 
ticle density of the different targets is between 1.36 and 
1.88 •IQ-'^* • cm~^ which corresponds to an areal density 
between (318 — 440) /ig/cm^. To study the homogeneity 
of the targets, the areal particles densities of several tar- 
gets were determined at many points distributed over the 
target area. An example is displayed in Fig. [31 A second 
RBS measurement after the activation experiment has 
shown that within the uncertainties, no target material 
was lost during the activation runs. 

The metallic Pr targets oxidize completely within a 
few days. Because the oxide is not bound on the backing 
anymore, it is mandatory to avoid oxidation before the 
activation and to fix the target material at its position 
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FIG. 4. Activation setup at the cyclotron of the PTB Braun- 
schweig. The negatively charged diaphragm {Us ~ —300 
V) suppresses secondary electrons to ensure that the whole 
charge is measured on target. To protect the target against 
decomposition it is water cooled from the outside. 



if oxidation occurs. Therefore, the targets were stored 
in vacuum before the irradiations. After the irradiations 
they were put in a paraffin envelope. 



B. Experimental installation at PTB 

Doubly charged He^+ ions were extracted from the cy- 
clotron to irradiate the ^^^Pr targets in an activation 
chamber, which is designed as a Faraday cup. The charge 
deposited on the target was recorded in time steps of 
60 seconds by a current integrator for later correction 
of beam-current fluctuations. The beam-current inte- 
gration at this setup is well established and an uncer- 
tainty of 1% has been taken into account. Secondary 
electrons were suppressed by a negatively charged di- 
aphragm {Us — —300 V) at the entrance of the activa- 
tion chamber. A homogenous illumination of the target 
was achieved by wobbling the a beam. For each energy 
the wobbling was optimated by inserting a quartz win- 
dow at the target position and checking its homogeneous 
irradiation. The beam spot was about ten mm in diam- 
eter. The backings were water cooled from the outside. 
In Fig. |4] a schematic of the target chamber at PTB is 
shown. The energy of the a-beam was defined within an 
uncertainty of ±25 keV by means of the field calibration 
of two analyzing magnets as well as by a time-of-fiight 
measurement of the particle velocity [4g| . 



C. Irradiation and 7 counting 

In total ten ^''ipr samples were irradiated at eight dif- 
ferent a energies between 11 MeV and 15 MeV to produce 
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the unstable reaction product -'^^^Pm. In order to test 
the thermodynamical stabihty of the targets, two targets 
were irradiated with different beam currents between 0.08 
and 3.5 /lA at a energies of 12.6 MeV and 15 MeV. These 
measurements yielded current-independent cross sections 
within the uncertainties, excluding decomposition of tar- 
get material. The Q value of the ^^^Pr(Q;,n)^'*'*Pm re- 
action is Q = (-10246.19 ± 2.70) keV and the reaction 
threshold amount to (10537.23 ± 2.78) keV [49|. The du- 
ration of the activation runs was varied between one hour 
and 17 hours. The average beam current was 3 /xA. 

The produced activity of all targets was measured off 
beam at the Institut fiir Kernphysik in Cologne using 
a clover-type HPGe detector with a relative efficiency 
of 120% at = 1332.5 keV, compared to a 3 x 3 inch 
cylindrical Nal detector. The activity was determined by 
performing spectroscopy of the 7 rays emitted after the 
electron capture of ^^*Pm with a low-background count- 
ing setup. The clover-type HPGe detector is composed of 
four crystals, which provides the possibility of measuring 
coincidences between the crystals. Since the (Q;,n) cross 
section at i?a= 11 MeV is in the range of 10 /ib only 
and the reaction product ^''^Pm decays with a long half- 
life of ri/2 = (363 ± 14) days H^, the count rate in the 
singles spectra is below or close to the sensitivity limit. 
Therefore, the requirement of coincidences is mandatory 
to enhance the peak-to-background ratio. The duration 
for the counting lasted between one and 21 days. In 
the following, this detector is referred to as the Cologne 
clover. 

In order to exclude systematic errors in the 77 co- 
incidence method, singles spectra of five targets were 
recorded additionally at PTE with a coaxial HPGe de- 
tector with a relative efficiency of 70% (PTE 70% detec- 
tor). The Cologne clover as well as the PTE 70% detector 
were used in a close geometry between the target and the 
detector end cap. Moreover, the decay of the reaction 
product involves three cascading 7 transitions. There- 
fore, coincidence summing effects, which occur when two 
or more 7 rays are recorded in one crystal within the 
resolving time of this crystal, have to be taken into ac- 
count when using a close geometry [sO]. Therefore, the 
efficiencies of these detectors were corrected for coinci- 
dence summing effects (see Sec. IIV A"|) . 

To verify this correction, a further measurement was 
performed at the Institut fiir Kernphysik in Cologne with 
a large distance between end cap and target by using 
a coaxial HPGe detector with a relative efficiency of 
55%. Due to this large distance, summing effects play 
a marginal role. This detector is named Cologne detec- 
tor in the following. 

Furthermore, it was possible to perform the spec- 
troscopy of the target which was irradiated with an en- 
ergy oi Ea — 11.4 MeV at the underground laboratory 
for dosimetry and spectrometry (UDO) at the Asse salt 
mine near Eraunschweig. PTE operates various low- 
background 7-ray spectrometry systems at this under- 
ground laboratory at a depth of 1200 m water equivalent. 
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FIG. 5. Comparison between the experimentally determined 
efficiencies without summing effects and the simulation for the 
Cologne clover. It can be seen that the Geant4 simulation 
describes the experimental efficiencies very well if summing 
effects are negligible. Only the energy region of interest is 
shown here and the energy region where the transitions in 



Nd occur is marked by the shaded area. 



A coaxial ultralow-background HPGe detector with a rel- 
ative efficiency of 90% (ULE) was used for this measure- 
ment. Detailed information about the ULE detector and 
the UDO laboratory can be found in Ref. [5l|. This mea- 
surement was analyzed completely independently from 
the other measurements mentioned above. Two of the 
targets were measured in addition at the PTE with a 
coaxial HPGe detector with a relative efficiency of 50% 
(PTE 50% detector) and were analyzed with the same 
routine as the one which was counted at the ULE detec- 
tor. 

Within the uncertainties the measurements with dif- 
ferent detectors gave consistent results. 



IV. DATA ANALYSIS 

The reaction product ^'*'^Pm decays via electron cap- 
ture with a half-hfe of T1/2 = (363 ± 14) days to "^Nd 
[igj, see Fig. [S] Subsequent to the decay, the reaction 
product is in an excited state, which deexcites via the 
emission of 7 rays. The three strongest 7 transitions used 
for the analysis proceed in a cascade; their properties are 
summarized in Table ID 

The counts in the full-energy peak of a certain transi- 
tion Y{E^) are connected with the total number of decays 
of the reaction product AA'^ during the counting time by 
the following expression: 



Y{E^)^^[E^y{E^) 



^LIVE 
iREAL 



(3) 
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TABLE I. Decay data 


of the reaction product ^''^Prn. 


Only 


transitions which were 


used for the data analysis are 


listed. 


The decay parameters 


are taken from [4^. 




1 keV 




Mult. 


476.78 ± 0.03 


0.4378 ± 0.0199 


E2 


618.01 ± 0.03 


0.9850 ± 0.0298 


E2 


696.49 ± 0.03 


0.9949 ± 0.0002 
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FIG. 6. Simplified decay scheme of Pm. Only transitions 
which were used for the data analysis are shown. Spins, par- 
ities, energies and population probabilities le are taken from 

Si. 



The full-energy efficiency is denoted as e{Ej) and the ab- 
solute 7 intensity per decay of the mother nuclide ^^"^Pm 
is indicated by I^(E^). The correction |iiv^ takes the 
dead time or the measurement into account, which was 
smaller than 1 % for all targets. 

From the decay law another relation for the total num- 
ber of decays of ^*'*Pm during the counting time can be 
derived: 



AA^ = iVact [1 



-AAi„ 



(4) 



The total number of ^'**Pm nuclei at the end of the acti- 
vation is A'act. The counting time is denoted by Atmoas, 
whereas the time between the end of the activation and 
the beginning of the counting is denoted as At^ait- The 
quantity A is the decay constant. 

During the activation A^prod ^*^Pm nuclei are produced 
in total but, during the activation, a part of these nuclei 
decay so that, as mentioned before, A^act reaction prod- 
ucts are left at the end of the activation: 



A^act — yactA^prod 



(5) 



For a constant beam current during the activation the 
factor /act is obtained by 



/a 



l-e 



-AAta 



AiactA 



(6) 



The activation duration is denoted as Aiact- The beam 
current was recorded in time intervals of 60 s and an 
appropriate correction for small fluctuations was applied. 



A. Determination of detector efficiencies 

For the calculation of the reaction cross section the 
absolute full-energy efficiencies of the detectors have to 
be known. In total, five HPGe detectors were used for the 
spectroscopy as described in Sec. IIII CI The efficiencies 
of all detectors were determined and the procedure is 
presented in the following. 

Because of the low activity of the targets, the distance 
between the target and the end cap of the Cologne clover 
was only 5.9 mm. As mentioned before, the decay of 
^^^Pm involves three cascading 7 transitions, so that co- 
incidence summing effects have to be taken into account. 
This holds for any multiline calibration source with cas- 
cading 7 transitions. 

The absolute detector efficiency was measured in far 
geometry (10 cm distance between end cap of the detector 
and calibration source), where coincidence summing is 
negligible. In total, six calibration sources (^°Co,^'^''Cs, 
152Eu, "3Ba, "Co, and 226Ra) were used. 

In a second step, the efficiency was measured for all 
six calibration sources in close geometry as well. For 
the 7 transition following the decay of ^'^^Cs a conver- 
sion factor was calculated between the efficiency in far 
and in close geometries. This is a pure geometric fac- 
tor, because no coincidence summation occurs. Simula- 
tions with Geant4 [s^I confirm that the assumption of 
an energy-independent conversion factor is valid in an 
energy range between 200 keV and 2000 keV. Thus, the 
experimental efficiency curve determined in far geometry 
is shifted by this factor to yield the efficiencies in the 
close geometry without any summing effects. 

Monte Carlo simulations with Geant4 were performed 
for this close geometry and were compared to the exper- 
imental efficiencies without summing effects. The sim- 
ulation describes the experimental efficiencies very well. 
This is shown in Fig. [S] 

In a third step, we simulated the extended geometry of 
our targets but, within the uncertainties, the simulation 
delivers the same results for a point source and our target 
geometry. 

In the last step, the 7 cascade of ^'*'*Nd was imple- 
mented into the simulation to take summing effects into 
account. The simulation does not include 77 angular 
correlation between the three 7 rays. Calculations have 
shown that the effect of this 77 angular correlation on 
the efficiency is smaller than 1%. 

The counting with the PTB 70% detector was per- 
formed in the close geometry as well. Therefore, the same 
procedure as for the Cologne clover was used for the de- 
termination of the absolute full-energy efficiency of the 
PTB 70% detector. In this case only three calibration 
sources (®°Co, ^^^jt;^^ and ^'^^Cs) were used. 
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FIG. 7. Sensitivity of the 77 coincidence method. In the 
upper panel (a) the singles spectrum for a target irradiated 
with 11 MeV a particles is shown, whereas the lower panel 
(b) shows the corresponding coincidence spectrum for a gate 
triggered by 696 keV photons. 



The efficiencies of the ULB detector and the PTB 50% 
detector were simulated with Gespecor [SJ]. This 
program uses recent data from the data center DDEF 
[55| . Summing corrections were also calculated with spe- 
cial routines of this code. For the rest analysis of the 
ULB spectra the program "Genie 2000 Gamma Analy- 
sis Software" by CANBERRA was applied, whereas the 
data of the PTB 50% detector were analyzed by using 
the program Interwinner by Ortec. 

Finally, a large distance of 9.5 cm between end cap 
and target was used with the Cologne detector. There- 
fore, summing effects play a marginal role and simula- 
tions with Geant4 were unnecessary. An interpolation 
between the measured efficiencies using two calibration 
sources ("'^^^Cs and -'^^^Eu) was sufficient to obtain the 
full-energy efficiencies at energies of the ^■*"'Pm emission 
lines. 



B. 77 coincidence method 

As mentioned above, the activity of the targets irradi- 
ated with a particles at 1 1 MeV is below or of the order 
of the background activity. Therefore, even the strongest 
transitions in ^'^''Nd are only weakly visible in the singles 
spectra of the Cologne clover. A 77 coincidence method 
was applied to suppress the background. 

All coincidences between two crystals in the list mode 
data were sorted into a coincidence matrix. This symmet- 
ric matrix contains all coincidence events of any crystal 
pairs. An example for the analyzing power of the coinci- 
dence technique is depicted in Fig. [T] 

In Fig. ^a) the singles spectrum of the target irra- 
diated with 11 MeV a particles is shown, whereas Fig. 
[Tljb) shows the corresponding coincidence spectrum after 



applying a gate triggered by 696 keV photons. Both spec- 
tra show the relevant energy range, where all three transi- 
tions in ^''^Nd are located. The energies of the transitions 
are marked in gray. In the singles spectrum the 7 lines are 
hardly visible and superimposed on a huge background. 
In the coincidence spectrum the transitions are clearly 
visible and the peak-to-background ratio is improved dra- 
matically compared to the singles spectrum. The broad 
lines in the coincidence spectrum result from Compton 
scattering of strong background lines from ^^*Bi, when 
696 keV of the total energy is deposited in one of the 
crystals and an energy of — -Etotai — 696 keV is de- 
posited in another one [56| . 

To determine the total cross section, the absolute full- 
energy efficiency for the coincidence spectrum is needed 
as well. Five targets were used, two of them were bom- 
barded at 15 MeV, the others were bombarded at 14.4 
MeV, 13.8 MeV, and 13.2 MeV, respectively. Conversion 
factors 77(£'gato, E) between the singles spectrum and the 
coincidence spectrum were then calculated by 



ViEgeitoyE) = 



-^single (^^) 
A'coin(-Egatc, E) 



(7) 



The yield in the singles spectrum at an energy E is de- 
noted -/Vsingio(£') and the yield in the coincidence spec- 
trum at an energy E gated by photons of an energy E'gate 
is Ncoin{Egate, E). For each of the three transitions two 
conversion factors were calculated depending on which 
coincidence pair is considered. In total, six conversion 
factors were obtained by averaging r]{Egs,tc, E) over the 
five targets. These conversion factors, which are inde- 
pendent of the target activity, were used to calculate the 
corresponding yields in the singles spectrum A'singic(£') 
for all targets. These yields are used to derive the cross 
sections for the Cologne clover as described in Sec. IIVI 



V. RESULTS AND DISCUSSION 

The experimentally determined cross sections obtained 
with each detector used in this experiment are shown in 
Table HT] The values for the different detectors are in 
excellent agreement. This proves that the 77 coincidence 
method can be used to determine absolute values for cross 
sections in activation experiments. 

The energy Ea of the a particles was obtained by cor- 
recting the adjusted primary energy Eq of the a-particle 
beam with the energy loss AE inside the target material 



Ea — 



E„ + {Eq - AE) 



(8) 



This energy determination, which is appropriate if the 
energy loss and the cross sections change only slightly 
over the target thickness, is independent from a specific 
cross section prediction. In this experiment the uncer- 
tainties of the cross section are larger than the changes 
in the cross section over the target thickness. The widths 
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TABLE II. Summary of the experimental cross sections listed for each energy Ea of the a particles together with the detectors 
used. The areal particle density m of the targets is also indicated. 



keV 


m 
cm~^ 




Detectors Used 


(7 

mb 


11047 ± 28 


(1.32 ±0.05) 


.10^« 


Cologne clover 


0.008 ±0.001 


11305 ± 28 


(1.42 ±0.06) 


lOi* 


Cologne clover 
ULB 

PTE 50% detector 


0.015 ±0.001 
0.015 ±0.001 
0.014 ± 0.001 


11899 ± 29 


(1.71 ±0.07) 


•10^** 


Cologne clover 


0.07 ±0.01 


12528 ± 29 


(1.80 ±0.07) 


.10i« 


Cologne clover 


0.31 ±0.02 


12549 ± 28 


(1.36 ± 0.05) 


• 10 *^ 


Cologne clover 
PTE 70% detector 


0.35 ± 0.02 
0.35 ± 0.03 


13098 ± 29 


(1.88 ±0.08) 


•10^** 


Cologne clover 
PTE 70% detector 


1.20 ±0.07 
1.20 ±0.10 


13736 ± 28 


(1.56 ±0.06) 


.10^« 


Cologne clover 
PTE 70% detector 


3.58 ± 0.24 
3.56 ±0.31 


14426 ± 29 


(1.61 ±0.06) 


•10^** 


Cologne clover 
PTE 70% detector 
r ± £> ou/0 detector 
Cologne detector 


11.02 ±0.68 
10.95 ±0.91 

Q QB -1- n (KC\ 

y.yo It u.ou 
11.40 ±0.88 


14956 ± 28 


(1.72 ±0.07) 


10^** 


Cologne clover 
PTE 70% detector 


23.30 ± 1.37 
23.69 ± 1.96 


14955 ± 28 


(1.65 ±0.07) 


.10i« 


Cologne clover 
PTE 70% detector 


26.29 ± 1.68 
24.07 ± 2.00 



of the energies is determined by the straggling in the tar- 
get and the uncertainty of the beam energy of the ac- 
celerator. They were added according to Gaussian error 
propagation. As mentioned in Sec. IIII Bl the energy un- 
certainty of the cyclotron at the PTB is ± 25 keV. The 
energy loss of the a particles in the target itself is ob- 
tained by a Geant4 simulation and varied between 62 
keV and 85 keV. The simulation yields the distribution 
of the a-particle beam after traveling through the target 
material. The maximum of this distribution is Eq — A_B 
and the halfwidth of this distribution at 1 /e of the maxi- 
mum is the straggling in the target which varied between 
13 keV and 15 keV. 

Table IIIII gives the astrophysical S factors obtained 
with weighted averages of the cross sections from the dif- 
ferent detectors shown in Table HIl Note that the energies 
in Table IIIII arc given as center-of-mass (cm.) energies 
(as in Figs. [TJ [21 El andE]) to facilitate comparison with 
theoretical calculations. Since the energy enters the cal- 
culation of the S factor, the inclusion of the uncertainties 
on cross sections and energies is not straightforward. A 
pair of S factors is given for each cm. energy Ec.m. in 
Table IIIII evaluated at the lower and upper limit of the 



energy range defined by the uncertainties i?c.m. i Ai?. 
The error bars on a energy and cross section translate 
into an error trapezoid for the S factor, with its four cor- 
ners given by the upper and lower limit of the S factor 
pair at each cm. energy. This error trapezoid is also 
shown in Figs. [5] and [HI 

As discussed in Sec. [Ill the cross section of the (a,n) 
reaction is almost exclusively sensitive to a change in 
the a width (or equivalently the total a Hauser-Feshbach 
transmission coefficient). Entering the calculation of the 
a width are the transmission coefficients for the individ- 
ual a emissions to the ground and excited states of i^ipj. 
which are computed by solving radial Schrodinger equa- 
tions using an a+nucleus potential and the appropriate 
quantum numbers of the involved states (37j . Therefore, 
the excited states in i^^Pr also have to be known or a nu- 
clear level density where the discrete states are unknown. 
In the laboratory cross sections and S factors discussed 
here, the a widths are dominated by transitions to low- 
lying excited states because transitions to levels at higher 
excitation energies have lower relative a energy and are 
strongly suppressed by the Coulomb barrier. The proper- 
ties of the relevant states are mostly known and therefore 
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FIG. 8. Astrophysical S factors for reaction ^''^Pr(Q,n)^''''Pm 
as function of cm. energy. The experimental S factors were 
computed from averaged cross section values from all detec- 
tors. These are compared to calculations with the codes 
SMARAGD [13, 1131 and TALYS [H, [H using a+nucleus 
optical- model potentials from ^ (McFV_p3] (AA), [U,!!! 
(RF), and the standard potential [H, iMpused in TALYS 
(STD). 
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FIG. 9. Astrophysical S factors for reaction ^''^Pr(Q,n)^''''Pm 
as function of cm. energy. The experimental S factors were 
computed from averaged cross section values from all de- 
tectors. These are compared to calculations with the codes 
SMARAGD [STII using the a-l-nucleus optical-model po- 
tentials from [231 (McF) and a potential with an energy- 
dependent depth of the imaginary part (fit) . 



TABLE III. Astrophysical S factors as a function of center-of- 
mass energy, calculated from the cross sections averaged over 
all detectors and evaluated at the lower and upper edges of 
the energy uncertainties Es = -Ec.m. ± Ai? from the averaged 
uncertainties. The uncertainties AS on the 5* factors at each 
energy Ec.m. ± Ai5 are derived including the averaged uncer- 
tainties of the cross section. The uncertainties are quoted for 
a one a confidence inteval (coverage factor k=l). 



E 


AE 


E^ 


S factor 


AS 


MeV 


MeV 


MeV 


10^^ MeVb 


10^^ MeVb 


10.742 


n 0981 8 












10.714 


3.0349 


3701 






10.770 


2.5375 


30945 


10.993 


n 09898 












10.965 


2.4108 


91 91 6 






11.021 


9 0989 


18438 

U . -L 04:00 


11.570 


00989 












11.541 


3.1276 


0.21116 






11.599 


2.6545 


1 7999 


12.202 


n 0979fi 












12.174 




n DQl ^7 






12.230 


1.6213 


0791 


12.182 


n 0991 fi 












12.153 


1.7686 


10898 






12.211 


1.5176 


09991 




u.uzoyu 












12.707 


1.6745 


0.10099 






12.765 


1.4531 


0.08764 


13.357 


0.02841 












13.329 


1.1414 


0.07502 






13.385 


1.0025 


0.06590 


14.028 


0.02864 












13.999 


0.78233 


0.04914 






14.057 


0.6929 


0.04352 


14.544 


0.02803 












14.516 


0.5853 


0.03425 






14.572 


0.52304 


0.03061 


14.542 


0.02841 












14.514 


0.63701 


0.04211 






14.570 


0.56837 


0.03757 



the only uncertainty arises from the a+nucleus optical- 
model potential (OMP). 

Figure [8] shows a comparison between the experimen- 
tal S factors from Table IIIII and calculations performed 
with the statistical model codes SMARAGD [sj [131 and 
TALYS [H [H using different a-|-nucleus OMPs. Both 
codes were used with their default settings; only the 
a+nucleus OMP was varied. This figure can be directly 
compared with Fig. [5] to judge the impact of the OMP. 
Only below 11 MeV - close to the (a,n) threshold - may 
variations in the 7 and neutron widths also influence the 
prediction to a small extent, as well as threshold effects 
like width fluctuations. Interestingly, none of the used 
OMPs can reproduce the energy dependence of the data, 
even above 11 MeV. 

The most widely used astrophysical reaction rates in- 
clude the OMP by [13] (denoted by 'McF' in Fig. [8]). 
This potential was obtained by fitting a four-parameter 
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Woods-Saxon potential on measured a-particle elastic 
scattering data at an energy of = 24.7 MeV on nuclei 
between oxygen and uranium. It was successful in repro- 
ducing a wide range of scattering and reaction data at 
higher energies but was found to systematically overpre- 
dict reaction data close to the astrophysically relevant 
energy region (see, e.g., [H, IH, IsTI - ISSi and references 
therein). Here, we find similar overpredictions from a 
factor of two in the region 11.5 < i?"^. < 12.5 MeV to a 
factor of four below 11 MeV. Only the data at the upper 
end of the measured energy range is reproduced satisfac- 
torily with this potential. This is interesting as it shows 
the transition from an energy region where the potential 
is appropriate to the region where it is not. 

By default, the TALYS code uses a simplification of 
the folding approach by [gl, [13] for the a OMP (denoted 
by 'STD' in Fig. IH). The resulting energy dependence of 
the S factors is slightly steeper than the one with p7| . 
It shows comparable discrepancy with the data at the 
lowest energies but reproduces the data at higher energy 
slightly better. It seems, however, as if the disagreement 
worsens again when continuing to even higher energies 
above the measured range. 

Motivated by the failing of the OMP by [13 , a modified 
potential was derived in Refs. [U [fi^l by simultaneously 
fitting low-energy reaction data in the A ~ 140 mass re- 
gion; more specifically, reaction data of ^''■^Nd(n,a)^^°Ce 
[62, i47Sni(n,a)i44Nd [65], and "4Sni(a,7)i48Gd [H. 
This OMP has the same imaginary part as [l^l but a 
shallower real part, resulting in generally lower cross sec- 
tions. It does not describe scattering data but that was 
never intended. The original notion was that the poten- 
tial should only be applied at astrophysically low energies 
(close to the Coulomb barrier) although an explicit en- 
ergy dependence of its parameters was not introduced 
due to lack of data. The comparison in Fig. |8] shows that 
indeed the higher-energy region is not reproduced well 
whereas the calculated S factors are closer to the exper- 
imental ones than those obtained with the other OMPs 
at lower energy. Interestingly, the energy dependence of 
the S factors calculated with this potential is very similar 
to the one obtained with the standard TALYS potential 
but the values are shifted in magnitude. 

The idea of [6l|, 111] to include reaction data into fits 
of the a OMP was recently picked up again by [60j . 
They fit mass- and energy-dependent potential param- 
eters to (a, 7), (a,n), and (a,xn) reaction cross sections 
closely around the Coulomb barrier for targets in the 
range 121 < A < 197 and also attempted to reproduce 
scattering data available at higher energy. This OMP in- 
cludes a volume and a surface term in the imaginary part 
(both mass- and energy- dependent) and has 48 parame- 
ters in total, which is to be c omp ared to the 7 parameters 
of the potentials by [13, HH, HJ. It has to be noted that 
the "'^"'^Pr(a,n) reaction data at higher energy by [iH!] were 
included in the fit. The general difficulty of extrapolating 
potential parameters to low energies is illustrated by the 
comparison of the S factors obtained with this OMP to 



TABLE IV. Parameters for the energy-dependent Q-|-nucleus 
optical-model potential derived from the data in this work. 





Real part 






Volume imaginary part 




'■c 


1.20 


fm 








V 

rn. 


= 185.00 
1.40 


MeV 
fm 


W = 


25 


MeV 
fm 


{l+cxp[(0.9Bc--Bc.m.)/2]} 

1.40 




0.52 


fm 


av = 


0.52 


fm 



our data in Fig. M (the OMP is marked by 'AA' in this 
figure). The three data points at the highest measured 
energies are reproduced perfectly but the S factors drop 
much too rapidly with decreasing energy below about 13 
MeV. 

As previously mentioned, the data in this work seem 
to cover regions where the standard potentials work ac- 
ceptably well and such where they fail. It can be at- 
tempted to bridge the two regimes by introducing a sim- 
ple ener gy- d ependence connecting the potentials by [13] 
and (6l1. |62| , following up on the original idea. The ex- 
ample of the OMP by [60], however, shows that one has 
to be cautious and that such extrapolations may not be 
generally valid. Therefore, we tried to limit the number 
of parameters as much as possible. They are listed in 
Table IIVI It was chosen to keep the depth of the OMP 
real part V and the geometry of real and imaginary part 
of [23] . Only the depth of the volume imaginary part W 
was made energy dependent. It has to approach the value 
of ^3 ^-t high energy but has to be shallower at energies 
below the Coulomb barrier energy Ec- To achieve this, 
similarly to [35|, a Fermi-type function was used: 



I _|_ e(0-9-Ec-Bc.„.)/aE ■ 

The relevance of Q.9i?c as the energy below which the 
potential parameters have to be strongly modified was 
pointed out in Ref. [60^]. The value = 2 MeV for the 
"diffuseness" of the Fermi-type function was varied but 
it was found that best agreement with the experimental 
data was achieved when using the same value as in Ref. 
[3^. An often overlooked fact is that the choice of the 
Coulomb radius parameter rc may have a strong impact. 
While the S factors obtained with the OMPs by [27| and 
[61I [6^ are quite insensitive to the chosen rc, the re- 
sults obtained with the energy-dependent OMPs by [11] 
and ^] are strongly affected (see also [13] for a further 
discussion). Also, the results with our energy-dependent 
OMP are sensitive to rc and we set rc = 1.2 fm for a 
best fit, which is quite similar to the 1.25 fm used in Ref. 

The S factors obtained with the fitted potential are 
shown in Fig. [9] along with the experimental values. 
For comparison, the values obtained with the energy- 
independent OMP by [13] are also included again. Al- 
though not all of the detailed features seen in the data 
are reproduced, the overall energy-dependence of the S 
factors is well described. (It is mentioned in passing that 
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the two data at the lowest energies should be given lower 
weight in such a comparison regarding the a OMP be- 
cause the a width is not the only uncertainty there.) 

As explained in Sees. |T] and [III the relevant quantities 
in the 7 process are the (7,0;) reaction rates [which, in 
turn, are calculated from (0,7) rates]. Therefore, it is in- 
teresting to examine the impact of our fitted OMP on a 
captures. A complete comparison as, for example, is per- 
formed in Ref. [63], is beyond the scope of this paper, but 
two cases can be discussed. Extending our OMP into the 
astrophysically relevant energy window for the reaction 
^'*^Pr(Q;,7)^*^Pm (see Fig. [T|) and calculating the stellar 
reaction rate for this reaction, it is found that the rate 
would be lower by three orders of magnitude compared 
to the one obtained with the potential of 27] . This is due 
to the very different energy dependence of the potential 
and the resulting cross sections. For comparison., the 
rate resulting from application of the OMP by [6l|, [g^ is 
lower by only 10% to 20% in this case. The reduction in 
the ^^^Pm(7,Q;) rate is the same as for the capture rate 
and thus it would become much slower than the com- 
peting ^'*^Pm(7,n). This removes the ^'*^Pm(7,Q;) deflec- 
tion point previously present in the 7-process path [66| . 
Nevertheless, astrophysically this is of minor importance 
because there are no stable seed nuclei of Pm to be pho- 
todisintegrated and it could only change the processing 
of downflows from higher masses, which are small. 

More relevant is the reaction ^''*Sm(Q;,7)^^*Gd, be- 
cause it determines the ^''^^Sm/^^'^Sm ratio in core- 
collapse supernova ejecta which can be determined via 
Nd/Sm ratios measured in meteoritic material [8], js^ [67|- 
[6|. A calculation of the S factors of ^^'^Sm{a,jy'^^Gd 
with our potential showed that it yields values compara- 
ble at low energies to the fit potential in Ref. [s^] (see, 
however, Ref. |37| for a cautionary discussion of further 
uncertainties in that OMP). Due to the different energy 
dependence, however, it yields different rates in the rele- 
vant temperature range 2 < T < 3 GK. Compared with 
the calculation shown in Ref. 35| the rate is higher by 
50% at 3 GK, which coincides with the rate given there 
at 2.5 GK, and is lower by a factor of 3 at 2 GK. Since the 
ejected ^*^Sm/^^^Sm ratio is determined by photodisin- 
tegration at the lower end of the temperature range, this 
would result in an even larger ^^^Sm/^^'*Sm ratio as re- 
ported in Ref. [HI . It would be very difficult to reconcile 
such a ratio with the ones found in meteorites. It is im- 
portant to note that the energy dependence of our OMP 
parameters may be more complicated than assumed and 
the use of the potential at lower energies than the ones 
measured in this work may not be warranted. 



cross sections. This method improved in this case the 
peak-to-background ratio in the coincidence spectra by 
a factor of more than 50 compared to the singles spec- 
tra, which minimize the sensitivity limit dramatically. 
Therefore, the 77 coincidence method allows measuring 
at lower energies close to the astrophysical relevant en- 
ergy range. A comparison with 4 further detectors proves 
that the 77 coincidence method is an excellent tool to in- 
vestigate cross sections down to the microbarn range. 

It was possible to measure close to the astrophysically 
relevant energy range, allowing us to study the validity 
of a-fnucleus optical-model potentials at low energies. 
A comparison of the experimental results with Hauser- 
Feshbach statistical model calculations confirms that the 
theoretical description of sub-Coulomb a transitions re- 
mains problematic. Especially the reproduction of the 
energy dependence of the excitation function by the theo- 
retical models is worse than in previously studied compa- 
rable cases. This may be connected to the closed neutron 
shell of ^^^Pr because strong discrepancies between ex- 
periment and predictions were also found for the reaction 
"'^*^Sm(Q;,7)^*^Gd in a previous study. 

A local energy-dependent optical potential was derived 
to improve the description of the present data, thereby 
proving the possibility to do so by just employing a suit- 
able potential. It was further applied to calculate the 
astrophysical reaction rates for ^'*^Pr(a,7), "'^''''^Sm(Q!,7), 
and their inverses. The revised rates would remove 
a 7-process path deflection at ^""^Pm and increase the 
^'*^Sm/^'*^Sm ratio produced in core-collapse supernovae, 
respectively. The energy-dependence of the potential, 
however, is still unconstrained in the energy range actu- 
ally relevant for the astrophysical 7 process and deflnitive 
conclusions are premature. 

The results underline the difficulties to determine an 
a-l-nucleus OMP which is globally applicable not only 
for different target nuclei but also at astrophysically rel- 
evant, low interaction energies. This also underlines the 
importance of measuring as close as possible or resonable 
to the astrophysically required energies. In our case we 
decided not to measure at even lower energies because 
here the cross section is also sensitive to the neutron and 
7 width and a determination of the Q;-|-^''^^Pr OMP is 
challenging. In general, further investigations are neces- 
sary and global improvements cannot be expected before 
additional low-energy data for a wider range of targets 
are accumulated. 



ACKNOWLEDGMENTS 



VI. SUMMARY AND CONCLUSION 

The cross sections of the reaction ^^^Pr(a,n)^^^Pm 
were measured between 11 MeV and 15 MeV. The 77 co- 
incidence method was used within one clover-type HPGe 
detector in an activation experiment to determine total 



The authors acknowledge the help of the PTB ac- 
celerator staff participating in this experiment and M. 
Schmiedel from the activity department of the PTB. 
Moreover, we thank I. Dillmann and V. Derya for fruitful 
discussions during the preparation of the experiment and 
the data analysis. Many thanks to N. Warr for his helpful 
corrections. We thank P. Mohr for his valuable comments 



12 



on the manuscript. ME, AH, and AS are members of the Forschungsgemeinschaft under the contract ZI 510/5-1. 
Bonn-Cologne Graduate School of Physics and Astron- TR is supported by the European Commission within 
omy. This project has been supported by the Deutsche the FP7 ENSAR/THEXO project. 



E. M. Burbidge, G. R. Burbidge, W. A. Fowler, and F. 
Hoyle, Rev. Mod. Phys. 29, 547 (1957). 

A. G. W. Cameron, Pub. Astron. Soc. Pac. 69, 201 
(1957). 

J. J. Cowan, F.-K. Thielemann, and J. W. Truran, Phys. 
Rep. 208, 267 (1991). 

G. Wallerstein et al, Rev. Mod. Phys. 69, 995 (1997). 
M. Arnould, S. Goriely, and K. Takahashi, Phys. Rep. 
450, 97 (2007). 

F. Kappeler, R. Gallino, S. Bisterzo, and W. Aoki, Rev. 
Mod. Phys. 83, 157 (2011). 

J. Audouze and J. W. Truran, Astrophys. J. 202, 204 
(1975). 

M. Arnould and S. Goriely Phys. Rep. 384, 1 (2003). 
M. Arnould, Astron. Astrophys. 46, 117 (1976). 
S. E. Woosley and W. M. Howard, Astrophys. J. Suppl. 
36, 285 (1978). 

M. Rayet et al, Astron. Astrophys. 298, 517 (1995). 

S. E. Woosley, D. H. Hartmann, R. D. Hoffman, and W. 

C. Haxton, Astrophys. J. 356, 272 (1990). 

C. Arlandini et al., Astrophys. J. 525, 886 (1999). 

Zs. Nemeth et al, Astrophys. J. 426, 357 (1994). 

T. Rauscher, A. Heger, R. D. Hoffman, and S. E. 

Woosley, Astrophys. J. 576, 323 (2002). 

S. E. Woosley and A. Heger, Phys. Rep. 442, 269 (2007). 

T. Rauscher, Proceedings of Science PoS(NIC XI)059 

(2011). 

C. Travaglio, F. K. Ropke, R. Gallino, and W. Hille- 
brandt, Astrophys. J. 739, 93 (2011). 
W. M. Howard, B. S. Meyer, and S. E. Woosley, Astro- 
phys. J. Lett. 373, L5 (1991). 

W. M. Howard and B. S. Meyer, in Proc. 2nd Int. Conf. 
Nuclei in the Cosmos, eds. F. Kappeler, K. Wisshak (lOP 
Publishing, Bristol 1993), p. 575. 

M. Kusakabe, N. Iwamoto, and K. Nomoto, Nucl. Phys. 
A758, 459 (2005). 

M. Kusakabe, N. Iwamoto, and K. Nomoto, Astrophys. 
J. 726, 25 (2011). 

W. Hauser and H. Feshbach, Phys. Rev. 87, 366 (1952). 
T. Rauscher and F.-K. Thielemann, At. Data Nucl. Data 
Tables 75, 1 (2000). 

T. Rauscher, Phys. Rev. C 73, 015804 (2006). 

W. Rapp et al, Astrophys. J. 653, 474 (2006). 

L. McFadden and G. R. Satchler, Nucl. Phys. 84, 177 

(1966). 

G. G. Kiss et al, Phys. Lett. B 695, 419 (2011). 
C. Yalgin et al, Phys. Rev. C 79, 065801 (2009). 

S. Harissopulos et al, Nucl. Phys. A 758, 505 (2005). 

Gy Gyiirky et al, Phys. Rev. C 74, 025805 (2006). 

Gy Gyiirky et al, J. Phys. G 37, 115201 (2010). 

W. Rapp et al, Phys. Rev. C 66, 015803 (2002). 

W. Rapp et al., Phys. Rev. C 78, 025804 (2008). 

E. Somorjai et al., Astron. Astrophys. 333, 1112 (1998). 

H. J. Brede et al., Nucl. Instr. Meth. 169, 349 (1980). 
T. Rauscher, Int. J. Mod. Phys. E., 20, 1071 (2011). 

R. A. Ward and W. A. Fowler, Astrophys. J. 238, 266 



(1980). 

[39] W. A. Fowler, G. E. Caughlan, and B. A. Zimmermann, 
Annu. Rev. Astron. Astrophys. 5, 525 (1967). 

[40] J. A. Holmes, S. E. Woosley, W. A. Fowler, and B. A. 
Zimmerman, At. Data Nucl. Data Tables 18, 305 (1976). 

[41] S. E. Woosley, W. A. Fowler, J. A. Holmes, and B. A. 
Zimmerman, At. Data Nucl. Data Tables 22, 371 (1978). 

[42] G. G. Kiss et al., Phys. Rev. Lett. 101, 191101 (2008). 

[43] T. Rauscher et al, Phys. Rev. C 80, 035801 (2009). 

[44] T. Rauscher, Phys. Rev. C 81, 045807 (2010). 

[45] M. Afzal Ansari, N. P. M. Sathik, D. Singh, and M. H. 
Rashid, J. Phys. Soc. Jap., 74, 1150 (2005). 

[46] G. G. Kiss et al, Phys. Rev. C 83, 065807 (2011). 

[47] Zs. Fiilop et al, Phys. Rev. C 64, 065805 (2001). 

[48] R. Bottger, private communication (2002). 

[49] NNDC Onhne Data Service, ENSDF database, 
http://www.nndc.bnl.gov/nndc/ensdf/ . 

[50] K. Debertin and R. G. Helmer, Gamma and X-ray Spec- 
trometry with Semiconductior Detectors (North Holland, 
New York, 1988).. 

[51] S. Neumaier, M. Woicik, H. Dombrowski, and D. Arnold, 
Appl. Radiat. Isot. 67, 726 (2009). 

[52] S. AgostineUi et al., Nucl. Instr. and Meth. A 506, 250 
(2003).. 

[53] O. Sima and D. Arnould, Appl. Radiat. Isot. 53, 51 
(2000). 

[54] O. Sima, D. Arnold, and C. Dovlete, J. Rad. Nucl. Chem. 
248, 359 (2001). 

[55] M. M. Be et al, NUCLEIDE, Table de radionucleide 
sur CD-ROM, Version 2-2004, CEA/BNM-LNHB, 91191, 
Gif-sur-Yvette, France. 

[56] S. Dabaneh et al, Nucl. Instr. Meth. A 517, 230 (2004). 

[57] T. Rauscher, code SMARAGD, version 0.8.3s (2011). 

[58] A. J. Koning, S. Hilaire, and M. C. Duijvestijn, in Pro- 
ceedings of the International Conference on Nuclear Data 
for Science and Technology - ND2004, Sep. 26 - Oct. 1, 
2004, Santa Fe, USA, MP, Vol. 769, edited by R. C. 
Haight, M. B. Chadwick, T. Kawano, and P. Talou 
(2005) p. 1154. 

[59] TALYS-1.2, User Manual, http:/7www.talys.eu 

[60] M. Avrigeanu and V. Avrigeanu, Phys. Rev. C 82, 
014606 (2010). 

[61] C. Frohlich, diploma thesis. University of Basel, Switzer- 
land (2002). 

[62] T. Rauscher, Nucl. Phys. A719 (2003) 73; erratum: 
Nucl. Phys. A725 (2003) 295. 

[63] S. Watanabe, Nucl. Phys. 8, 484 (1958).. 

[64] D. G. Madland, in Proceedings of a Specialists' Meet- 
ing on preequilibrium nucelar reactions, edited by 
B. Strohmaier (1988) p. 103. 

[65] Y. M. Gledenov et al, Phys. Rev. C 62, 042801 (2000). 

[66] T. Rauscher, Phys. Rev. C 73, 015804 (2006). 

[67] A. Prinzhofer, D. A. Papanastassiou, and G. A. Wasser- 
burg, Astrophys. J. Lett. 344, L81 (1989). 

[68] S. E. Woosley and W. M. Howard, Astrophys. J. Lett. 
354, L21 (1990). 



T. Rauscher, F.-K. Thielemann, and H. Oberhummer, 
Astrophys. J. Lett. 451, L37 (1995). 



